Is the nonlinear Meissner effect unobservable? 
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We examine the effects of nonlocal electrodynamics for a d-wave superconductor on the field depen- 
dence of the magnetic penetration depth. The linear field dependence predicted in the local limit, 
commonly known as the nonlinear Meissner effect, is instead found to be quadratic, 5 A ~ H 2 for 
fields below a crossover scale H*. This crossover is shown to be geometry dependent, and for most 
orientations of the screening currents, is of the same order as or greater than H c i, implying that 
the nonlinear Meissner effect can not be observed. For special orientations where the current flows 
along the nodal directions, however, the nonlinear Meissner effect may be recovered. 



PACS Numbers: 74.25.Nf, 74.20.Fg 

Measurements of the linear term in the temperature 
dependence of the electromagnetic penetration depth |l]] 
5X(T) = A(T) - A(0) - T played a pivotal role in the 
identification of the d-wave symmetry of the cuprate su- 
perconductors, ft was always clear, however, that they 
were sensitive only to the depletion of the superfluid frac- 
tion by the thermal excitations of quasiparticles with mo- 
menta near the order parameter nodes, and hence pro- 
vided information only on the energy gap on these regions 
of the Fermi surface. Furthermore, in a d-wave supercon- 
ductor with tetragonal crystal symmetry, the penetra- 
tion depth is independent of the angle 9 the supercurrent 
makes with the crystal axes since the linear response ten- 
sor is isotropic; the gap shape is thus not directly reflected 
in the angle dependence of the effect. 

Yip and Sauls argued, however, that information 
about the gap shape in an unconventional superconduc- 
tor with line nodes could be obtained by measuring the 
nonlinear (magnetic field-dependent) penetration depth 
SX(H) — X(H) — A(0), and predicted an unusual field 
dependence, 



SX(H) 



(1) 



where is a (weak) angle-dependent function, H is the 
applied magnetic field, and H = 3$o/( 7I ' 2 ^oAo) is of the 
order the thermodynamic critical field of the system (<£>o 
is the flux quantum, £o the coherence length and Ao the 
penetration depth). Later, it was shown || that im- 
purities alter this prediction, leading to an asymptotic 
SX ~ H 2 dependence, but only at concentrations such 
that the linear term in the H — > temperature depen- 
dence SX(T) ~ T is replaced by SX(T) ~ T 2 §-§. 

From an experimental standpoint, the search for con- 
firmation of this fundamental manifestation of d-wave 
symmetry has been controversial. In 1995, Maeda et 
al. reported the first observation of a linear H term on 
Bi 2 Sr 2 CaCu 2 8 {BSCCO -2212) films §, but the re- 
sults were questionable since they were obtained at rel- 
atively high temperatures, where the Yip-Sauls theory 
predicts a SX ~ H 2 . The clean d-wave nonlinear Meissner 
effect should be obtained, according to the simple theory, 



only for temperatures below a crossover T <C E non n n , 
with E non ii n = v s kp ~ (H/ Hq)Aq and Ao the gap maxi- 
mum, giving typically E non n n ~ IK for the cuprates, for 
fields close to H c \, the lower critical field. Measurement 
of a linear- term was also claimed at much higher tem- 
peratures for YBa 2 Cu?,0-i-s (YBCO) samples ||. It is 
important to recall that several extrinsic effects can give 
rise to spurious field dependences, most importantly con- 
tributions from trapped vortices which enter near sample 
edges. An extrinsic contribution 5X ~ >/~H to the effec- 
tive penetration depth can arise from coupling of the rf 
field used in resonant coil experiments to the pinned vor- 
tices 0. Such dependences were also found in measure- 
ments at much lower temperatures and on higher quality 
sing le crystals of YBCO by Carrington and Giannetta 
fllOf . Finally, a measurement to detect higher-harmonic 
nonlinear modes in the ac response predicted on the basis 
of Yip-Sauls theory by Zutic and Vails O], apparently 
less sensitive to the effects of trapped flux, failed to ob- 
serve the predicted signal | Q . 

We therefore take the point of view that the available 
experimental data do not support the existence of a non- 
linear Meissner effect in the best cuprate samples, and 
ask why. Our answer is that the effects of nonlocal elec- 
trodynamics, neglected in the work of Yip and Sauls, 
serve to cut off the (local) nonlinear Meissner effect as 
the field is lowered, replacing the linear- H response with 
a quadratic one below a crossover field H*. In most ex- 
perimental geometries, this crossover is of the order of the 
lower critical field H&, meaning the SX(H) ~ H behav- 
ior cannot generically be observed in the Meissner state. 
However, we show that nonlocal effects are negligible in 
geometries with currents parallel to a nodal direction. 

Nonlocal electrodynamics. The Yip-Sauls description 
of the nonlinear Meissner effect is performed in the lo- 
cal limit, i.e. the d-wave Cooper pairs, although ex- 
tended objects, are assumed to respond to the electro- 
magnetic field at a single point at the center of the pair. 
Such an approximation is normally justified in classic 
type-II superconductors, where Ao is much larger than 
£o- Although the cuprates are strongly type-II super- 
conductors, with small in-plane £o — 0(15A) and large 
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in-plane Ao — 0(1500 A), effects of nonlocal electrody- 
namics are to be seen in several special situations, as a 
result of the gap nodes, along which the "effective co- 
herence length" VF/nAk diverges |l|. Such effects were 
shown to significantly suppress the shielding supercur- 
rents in measurements of the ac Meissner effect [fl5| in 
unconventional superconductors with line nodes, with 
application to heavy fermion materials. In the context 
of the cuprates, it was shown by Zucarro et al. [jl6| that 
significant variations in the electromagnetic response of 
clean <i-wave superconductors were to be expected when 
nonlocal effects were included. But it was Kosztin and 
Leggett who first pointed out that the very signature 
of a superconductor with line nodes, the linear-T term 
in the dc penetration depth, was modified at the low- 
est temperatures |14 . The crossover temperature below 
which this occurs is E non i oc ~ k _1 Ao, where k = Xq/^q 
is the Ginzburg-Landau parameter. The characteristic 
crossover field defined by E nordoc = E nonHn is given by 
H* ~ K _1 i?o — H c ±. One might therefore expect that 
nonlocal effects "cut off" the nonlinear Meissner effect. 
This observation needs, however, to be supported by a 
full calculation, which we now sketch. 

General current response. The current response for a 
d-wave superconductor may be derived from a full cal- 
culation of the BCS free energy for fixed external field 
H corresponding to vector potential A. We have per- 
formed such a calculation (including both nonlocal and 
nonlinear effects analytically) under the assumptions of 
i) negligible spatial variations of the amplitude of the or- 
der parameter, Ak(r); and ii) slow spatial variations of 
the superfluid velocity v s . The details of the derivation 
will be given elsewhere jlTj, but the final result for the 
Fourier component of the current in an infinite system at 
constant v s is j(q) = — /C(q, v s , T)A(q), where 



E(q,v a ,r) 
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with Ao = \J mc 2 /Anne 2 the local penetration depth at 
T = in linear response theory, n the electron den- 
sity, m the electron mass, c the speed of light, kn the 
projection of k onto the surface, Zk = iuj n + v s • kp, 
El = £ + A|, £± = £k ±q /2, and £ k = ^ - fi (p 
the chemical potential). This is the standard linear re- 
sponse result with all quasiparticle Matsubara energies 
u> n modified by the semiclassical Doppler shift v s ■ k^ 
p8[ . For the real physical system in the Meissner state, 
a surface boundary is present and v s decays with the 
distance from the surface y. In the linear, local approx- 
imation, v s (y) — (eXoH/mc) exp(— y/\o) decays expo- 
nentially and v s (y = 0) is proportional to the exter- 
nal field, neither of which hold generally. But we may 
imagine the system to be subdivided into many layers, 
in each of which v s is roughly constant. Furthermore, 
we assume specular scattering surface boundary condi- 
tion on the quasiparticle wavefunctions. In this case, the 



response is given by Eq. (^|). The geometry which we 
study first is such that the magnetic field is along the c 
axis and thus v s lies in the ab plane; the normal of the 
boundary plane q is also in the ab plane, forming an angle 
8 with the b axis. We will discuss the situation where H 
is perpendicular to the c axis later. 

To a first approximation, we neglect the space de- 
pendence of v s and replace it by its value at the sur- 
face v s ~ v s (y = 0) = eX Q H/mc as given by the solu- 
tion to the linear, local electrodynamics problem. Now 
we separate out the T — local, linear response as 
/C(q,v s ,T) = c/(47rAg) + <5/C(q, v s ,T), and then define 
the nonlinear, nonlocal penetration depth to be 
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This expression has the virtue of reducing exactly to the 
nonlocal expression of Kosztin and Leggett jl3| if the v s - 
dependence is neglected, and (qualitatively) to the non- 
linear expression of Yip and Sauls |2j if the g-dependence 
is neglected. It is worthwhile reviewing the latter case. 
When T < D s i;f we obtain 



where 



5X ; 
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n(9| 3 . Thus, 
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We recall that Yip and Sauls g] and Xu et al. || de- 
fined penetration depth from its initial decay: A = 
—H(Q)/[dH(y)/dy]\ y= o which is equivalent to the defi- 
nition in Eq. (|J) only in the linear limit. Using their 
definition we find that the prefactor 3/2 in Eq. (^) is 
modified to 3. At this stage, we may improve our ap- 
proximation by restoring the spatially varying nature of 
v s . We follow Yip et al. to solve directly a nonlinear 
London equation 

V v s = -/C(0+,v s ,T)v s = ^(1- __„.), (6) 

under the boundary condition dv 8 /dy\ y= o — (e/mc)H. 
Using the relation dv s /dy = (e/mc)H(y) and their def- 
inition for the penetration depth again, we see that 
SX/Xq = So in the local limit we get qualita- 

tively the same nonlinear correction to the penetration 
depth as Yip and Sauls j|, comparing to Eq. (|l|), but 
with a larger prefactor. This deviation arises from the 
perturbative treatment of the response of the supercon- 
ductor to the A(q 7^ 0) modes in the present theory. 
The nonlinear effects obtained at T <C v s kp, due to the 
Doppler shift coming from the coupling of the A(q = 0) 
mode to the quasiparticles, are nonperturbative results. 
A discussion of this point can be found in Ref . n% . 
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Low-energy scaling properties. Explicit analytical re- 
sults for the full response in the approximation of con- 
stant v s may now be obtained by expanding SIC for low 
energies, v s kF,qVF,T <C Aq. A 2-parameter scaling re- 
lation may be derived, 



u z=±i 



T 1 T 



Fx(z 1 ,z 2 ) = -zi + [ln(e Z2 + 1) + ln(e" 



1)] 



dx[f(x - z 2 ) + f(x + z 2 )][l - (-) 2 ] 1 / 2 , (7) 
o z i 



with a = ypq/i^Vz), £ — v s kp/\/2, ugi — | cos 9+1 sin#|, 
Wei = \s'm8 — lcos9\ and f(x) = 1/(1 + e x ). We consider 
primarily the limit T <C a, e. It is straightforward to see 
that for z 2 ^s> 1 



F\(zi,z 2 ) 



Z2 

4 



z\ <C z 2 
zi{z\ + tt 2 )/(6z 2 ) Z 2 <Zl. 



(8) 



The scaling form ([?]) may now be inserted into (|^) to find 
the penetration depth. Its field dependence turns out 
to be determined by the ratio awgi/eugi, with a typical 
q set to 1/Ao in (|J) by the denominator. If we define 
two new dimensionless parameters hoi = (ugi/wgi)h, I = 
±1, with h = mv s \o — > (3/tt)kH/Hq, we get 5\ spcc = 



J2i=±i ^spec, where 



5\,. 



(I) 



AV2 



Ug,K 1 Il 



cgiin — 



^A 2 
for hoi > 1, (9) 



with cgn — 0.008w| ; and cgi 2 = 0.006wg ; /wg;, and 



SX, 



Ao 



8^2 



dgi\K 1 + dgi 2 K 



- X h 2 



^0 



for het^l, (10) 



with dgn = 0.5ug Z w 9 i, dg i2 = \.\u\Jwgi and dg a = 
0A7ul t /wgi. The singularity which leads to the linear 
term in Eq. (||) if nonlocal effects are simply neglected 
is seen to be cut off by the constant term in Eq. (10). 
This happens at hgi ~ 1. For 9 not too close to a nodal 
value, wgi and ugi are of order unity for both I = ±1, 
and the crossover field may be simply determined from 
h = 1, H* ~ k~ 1 Hq. The thermodynamical critical field 
i? c ~ ffo is the geometric mean \JH c \H c2 of the upper 
critical field H c2 — <I>o/(27r£ 2 ) and the lower critical field 
H c i = ^0/ (ZkAq). Thus the crossover field for such ge- 
ometries is of the same order as the lower critical field, 
H* ~ H cl . Since the Meissner state is unstable to the 
Abrikosov vortex state for fields above H c i , the linear- H 
term in 5X is asymptotically unobservable in these ge- 
ometries. In the Meissner state it is, in fact, replaced by 
a quadratic variation, as seen in ([To|). In Fig. 1, we have 
plotted the penetration depth calculated in Kosztin and 



Leggett's geometry 9 = 0, according to the exact scaling 
function given in (Q), to explicitly exhibit this effect. 

When 9 is near any node, however, ugi can be much 
larger than wgi and the nonlocal energy scale much 
smaller than the nonlinear one, leading to a very small 
crossover field H* ~ H cl woi/ugi <C H c \. We conclude 
that the nonlinear Meissner effect may be observable in 
such special geometries, e.g. a sample with a dominant 
(110) surface and H || c. 
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FIG. 1. Normalized penetration depth correction 8\/\q 
vs. normalized magnetic field H/Hq at T = for a 
sample with H parallel to a (010) surface (9 = 0) and 
Ginzburg-Landau parameter k = 10 2 , hence lower critical 
field H c i — W~ 2 Ho- Squares: full response; solid line: nonlo- 
cal effects neglected. Insert: T = 0.1T C . Circles: full response; 
solid line: nonlocal effects neglected. Dashed line: H* ~ H c \. 

There is another special situation in the quasi-2D 
cuprates in which the characteristic energy of nonlocal 
effects becomes very small and the previous discussion 
might not be expected to apply. If we assume that the 
quasiparticles are strictly confined to the ab plane and if 
the experimental configuration is dominated by a (001) 
surface with H || ab, then = q • Vf = 0. In 

this case, however, the magnetic field is not screened 
(A — > 00) Jl9|]. Whenever the quasiparticles are al- 
lowed to move along the c-direction we expect nonzero 

^nonioc Pol - This can be demonstrated directly for any 
model of coherent transport along the c-axis. The ma- 
terials of greatest current interest are the cuprate mate- 
rials YBCO and BSCCO - 2212. In the YBCO case, 
it is reasonable to treat the system as weakly-3D. The 
characteristic nonlocal energy is then vp c q ~ £o c A /Ao 
much smaller than in the (010) case since the c-axis co- 
herence length is £oc — 3A as opposed to the in-plane 
coherence length of £0 — 15A On the other hand, the 
lower critical field is also much smaller for this geometry, 



H 



(ab) 



$o/(27rA A 0c ). Using A 



Or 



0.5-0.8 xl0 4 A, we 



again find a large crossover field i/*( ab ) ~ (£,Oc/£,o)H* — 
0.5 — lH^ b \ making it still impossible from a practical 
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point of view to extract a linear- H term [[21) . 

The BSCCO system is so anisotropic that it is ques- 
tionable whether the above argument applies. A full 
treatment of this problem awaits a generally accepted 
theory of the (incoherent) c-axis transport in the normal 
state. 

In conclusion, with the exception of the claim by 
Maeda et al. all attempts to measure the nonlin- 

ear Meissner effect, a fundamental consequence of d-wave 
symmetry, have failed. As the order parameter symmetry 
is by now well-established in these materials, an expla- 
nation is required. To this end, we have considered the 
interplay of nonlinear and nonlocal effects in the elec- 
tromagnetic response of a d-wave superconductor in the 
Meissner regime. This situation differs from the simi- 
lar problem considered by Amin et al. p2| in the vortex 
phase due to the fact that the nonlinear and nonlocal en- 
ergies E non ii n and E non i oc are not independent in their 
case, as a consequence of the fluxoid quantization condi- 
tion. By contrast, in the Meissner state, while the nonlo- 
cal scale is set by intrinsic normal-state parameters, i.e. 
T c and k, the nonlinear energy E non i in ~ v s kp is effec- 
tively tunable by an external held. We have shown that 
the nonlocal terms in the response generically cut off the 
nodal singularity in the nonlinear response, eliminating 
the linear- H term in the penetration depth. 

Neither the current theory nor the Yip-Sauls result, 
both based on a consideration of the shielding currents 
in the Meissner phase, can explain the results of Maeda et 
al. Maeda's results were obtained on samples which did 
not clearly exhibit a 5\(T; H — > 0) ~ T behavior in any 
range, and therefore should not show nonlinear effects 
even according to the Yip - Sauls theory. Furthermore, 
they were performed at temperatures of 9K and above, 
and a linear- H behavior was reported up to 50K. As we 
show in Fig. 1, even for temperatures as low as 0.1T C , 
no linear term should be visible even if nonlocal effects 
are neglected. Therefore we must conclude that Maeda 
et al. observed an extrinsic effect, probably related to 
trapped vortices. Very recently, new high-resolution res- 
onant coil measurements @,^3| on very clean YBCO 
samples observed (in some geometries) linear H terms 
with magnitudes close to the Yip-Sauls prediction Jl0|,^3| , 
but simultaneously measured qualitatively inconsistent 
temperature dependence. These authors also observed 
effects attributed to trapped vortices |[(|, and we there- 
fore conclude that the linear term must be of extrinsic 
origin. Clearly, a test of the nonlinear effects which is 
independent of trapped flux is desirable. The harmonics 
of the nonlinear response predicted by Zutic and Vails in 
principle provide such a test, but to date, a high resolu- 
tion nonlinear transverse magnetization experiment de- 
signed to observe these resonances has failed to do so Q . 
This null result is consistent with the current theory, but 
by itself cannot confirm it due to uncertainties regarding 
the identification of H c \ 

Finally, in all experiments up to now the orientation of 
the supercurrent was along the crystal axes. Our analysis 



suggests the best possibility to observe the effect is in 
geometries with v s along the nodal directions. 
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